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Abstract. We semiclassically derive the leading off-diagonal correction to the 
spectral form factor of quantum systems with a chaotic classical counterpart. To 
this end we present a phase space generalization of a recent approach for uniformly 
hyperbolic systems 013- Our results coincide with corresponding random matrix 
predictions. Furthermore, we study the transition from the Gaussian orthogonal 
to the Gaussian unitary ensemble. 

PACS numbers: 03.65.Sq,05.45.Mt 

Advanced semiclassical methods have been very successful to significantly improve 
our understanding of complex, classically chaotic quantum systems 0]. This holds 
particularly true for observables which can be deduced from the (single-particle) Green 
function of the quantum system, such as the density of states, photo absorption, or 
orbital magnetism, to name a few. However, for quantities which are based on Green 
function products the situation is much more involved. These include linear response 
functions, e.g. for quantum transport, spectral correlation functions or, more generally, 
spectral statistics. A semiclassical treatment of such quantities is usually faced with 
the serious problem of evaluating multiple infinite sums over phase-carrying classical 
paths, which arise from the semiclassical representation of the Green functions in the 
limit h — > 0. A prominent example is the spectral two-point correlator or its Fourier 
transform the spectral form factor K(t). Random matrix assumptions leading to the 
prediction of a universal form for K(t) for classically chaotic quantum systems are 
supported by experimental and numerical data for a vast number of systems from 
different disciplines in physics [3]. Nevertheless, a corresponding analytical approach 
proving the conjectured 5 universality by explicitly including the underlying chaotic 
classical dynamics is still lacking. In this respect semiclassical techniques, bridging 
classical and quantum dynamics, appear to be natural tools but have to cope with the 
above mentioned problems when evaluating spectral correlations. For the form factor 
this involves the computation of (energy) averaged double sums over periodic orbits 
which has been addressed in several semiclassical approaches [H]. 

Recently, K(t) was semiclassically investigated for uniformly hyperbolic two- 
dimensional systems where the classical dynamics is governed by a single Lyapunov 
exponent By going beyond the usual diagonal approximation (K^^t) « 2r, 

describing the limit of spectral long-range correlations), the next to leading order 
contribution to K(t) was computed for systems with time-reversal symmetry. This 
was achieved by identifying off-diagonal pairs of correlated periodic orbits which are 
associated with each other via selfcrossings in configuration space. Based on this orbit 
class the random matrix theory (RMT) prediction for the form factor in the Gaussian 
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Orthogonal Ensemble (GOE), K(t) w2r- 2r 2 for r -> 0, could be derived PE], see 
also However, the important question remains whether the above result, as well 
as the RMT prediction, are specific for systems with uniformly hyperbolic dynamics, 
or whether they pertain for the much broader class of chaotic systems with different 
periodic orbits having different Lyapunov exponents. 

Here we present a generalization of the semiclassical approach outlined above to 
such non-uniformly hyperbolic systems in two dimensions and show that under rather 
general conditions the term ~ — 2r 2 in K{t) is indeed retained. To this end we develop 
a canonically invariant approach which is based on phase space arguments only. We 
identify as the relevant objects 'crossing regions' in phase space which can involve 
more than one selfcrossing in configuration space. We express the action differences 
of the considered orbit pairs in terms of local phase space properties, the directions 
of the stable and unstable manifolds, and present a method for counting the 'crossing 
regions'. This allows us eventually to determine their contribution to K(t). 

The semiclassical limit implies a large energy E compared to the mean level 
spacing l/d(E), The energy dependence of K(t) is smoothed out by an average over 
a classically small but quantum mechanically large energy window of size AE. In the 
considered limit AE can be chosen such that l/d(E) <C AE <C E. In this energy 
regime all classical actions, such as the action S spo of the shortest periodic orbit, are 
much larger than h, i.e. S spo /h 3> 1. Then one can employ Gutzwiller's trace formula 
|3] for the oscillating part of the density of states and evaluate the Fourier transform 
of the spectral two-point correlation function. This gives for the form factor a double 
sum over periodic orbits 

K(r) = ^yjA 1 AU^(i^)s(T-^±^.)\ , (1) 



AE 



with the time scaled according to t = T/Tn where Th = 2-KhdiE) is the Heisenberg 
time. We represent each periodic orbit 7 in terms of its phase space coordinates 
x = {QiP) where q and p are the position and momentum coordinates in two 
dimensions. In T 7 is the period of an orbit x~ ( and A 1 includes both, its weight and 
Maslov index. The action difference between two orbits x~ t and Xj is given by S-y^. 
The argument of the double sum is a rapidly oscillating function and hence the energy 
average suppresses most of the terms. To obtain a nonvanishing contribution to K(t) 
the action difference must be small, i.e. S 7 .j < h. Thus, for the time reversal case, 
the largest contribution is due to pairs of one path x 1 with itself or its time-reversed 
partner x 7 = (q 7 (T 7 — t), — p 7 (T 7 — t)) with vanishing action difference. Including 
only this type of pairs, known as the diagonal approximation [7j, reproduces the linear 
contribution K^(t) w 2r. 

In the recent approach beyond the diagonal approximation pairs of closely related 
periodic orbits with a small action difference have been constructed by analyzing self- 
crossings of the orbits in configuration space PE]. Here we proceed in a different way 
by generalizing the approach via self-crossings to 'crossing regions' in phase space: We 
will show that if a periodic orbit x 1 comes close to its time-reversed version £c 7 in such 
a 'crossing region' this can imply the existence of yet another periodic orbit a? 7 = x^ 
with a small action difference <S 7l7 between the two. The partner orbit x? follows the 
original a? 7 in a first segment 1Z and then the time-reversed path £C 7 during the rest 



of the time in the second segment C, see figure 1(a) It will turn out that not every 



'crossing region' implies the existence of a partner orbit and therefore a small number 
of non-relevant 'crossing regions' has to be excluded. 
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Figure 1. (a) Sketch of a correlated orbit pair in phase space (shown is a pro- 
jection of the four-dimensional space). The original periodic orbit x, the time 
reversed orbit x % , and the partner orbit x v are represented by the solid, dashed 
and dotted line, respectively. Due to time reversal symmetry each 'crossing region' 
appears twice. The Poincare surface of section (PSS) defined by the perpendicular 
coordinates (Ax- 1 -, Ap ) is indicated in the upper 'crossing region', (b) PSS with- 
in the 'crossing region' at the beginning of loop TZ. The path x defines the center 
of the coordinate system while the vector 8y points towards the time-reversed 
orbit x l . The position of the partner orbit x v is given by the vector SxR.i- 



To compute the contribution of the described orbit pairs (x 7 ,xP) to K(r) we 
rearrange Q into a sum over periodic orbits and another sum over all the partners 
with small action difference. This is based on the assumption that the dominant off- 
diagonal contribution to K{r), equation JJJ, is due to the systematic correlation of 
actions of the considered orbit pairs (x 7 ,xP) while other correlations are negligible. 
Then we sort the terms in the sums with respect to their action differences. Since 
the orbit length T 7 is proportional to Th for fixed r, one expects a large number 
of 'crossing regions' for each orbit. This allows us to replace the sum over action 
differences by an integral. The first off-diagonal contribution to K(t) then reads 



K (2) {t) = 4rRe 



dS 



/ dN s , 
\ dS 



(7,T) 




(2) 
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where A/g )7 is the number of relevant 'crossing regions' for a given periodic orbit x 1 
with an associated action difference smaller than S. In equation J2J, (• • -}( 7 t) denotes 
a weighted average over all orbits x 1 of given length T 7 = tTh- ft is defined as 

HA| 2 <5(T-T 7 ). (3) 



/ dN Sn 
\ dS 



(7,T) 1 7 



dS 



We proceed with the evaluation of © by first determining the geometry of the partner 
orbit xP. Then we show that the action difference between x 1 and x? is indeed small 
if the orbit x 1 and its time-reversed version x^ come close together in parts of the 
phase space. Finally we derive the averaged number (dNs n /dS) , T =T j of relevant 
'crossing regions' in the contributing regime S < h. 

Partner geometry. To construct the partner orbit we analyse the linearised 
equations of motion around x n in part TZ and around x 1 ^ in section C in the Poincare 
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surface of section (PSS), see figure^ an d show that a nontrivial solution representing 
the partner orbit x? exists under certain conditions. The distance between the 
original orbit x 1 and its time-reversed partner x^ in the PSS defined by the local 
transverse coordinates at the phase space position x = x^(t) is given by the vector 
Sy = (Ax- 1 , Ajr 1 -) ^U]. In general, the time evolution of a small deviation is determined 
by the starting point Xq of the original path in phase space and the initial deviation in 
the PSS Syo, e.g. Sy — 8y(t; xq, Syo)- Let us assume for the moment that the 'crossing 
region' under consideration is characterized by a small distance Sy. Decomposing 
8y = yuW u + y s w s in terms of the local unstable and stable manifolds uJ u . s !.) with 
the expansion coefficients 

Vu.s = ^ T „ z and Z = (4) 

the smallness of Sy is given if \y us \ <C 1. For simplicity we choose the relative 
orientation and lengths of the w UtS such that (Zw s ) T w u = S spo . It will turn out 
that \y u ,s\ <C 1 is the relevant regime for the evaluation of K(t) in the semiclassical 
limit. The distance between x 1 and the partner orbit x? at the beginning of the first 



loop 1Z is denoted by 8Sr,%, see figure 1(b) This vector lies in the PSS defined at 
the phase space position x = Xy(t) before the loop 1Z [this corresponds to the upper 



crossing region in figure |T(a)] . Having passed loop 1Z after time Tr this distance has 
changed to R 8xR_i with R being the stability matrix for loop 1Z. Before (and after) 
the other part of the orbit the difference between the time-reversed path x^ and the 
partner x? is denoted by Sx^ij (and L l Sx^i) where L % = FL~ 1 F is the stability 
matrix of the time-reversed loop C. The matrix F is defined as 

*-(; -,)■ <*> 

Solving the linearised equations of motion under the condition that the two parts of 
the partner orbit fit together in the 'crossing regions' yields the geometry of x? in 
terms of the distance Sy between the original orbit x 1 and its time-reversed xl f : 

Sx R<i =[l-L i R]~ 1 [l-L i F]6y. (6) 
The corresponding condition for Sx^ i for the lower 'crossing region' in figure 1(a) is 



found in a similar way This set of solutions Sx includes terms of order 0{y u ^ s ) 
and defines the partner orbit for a given small Sy representing a 'crossing region'. 

We now argue that a 'crossing region' does not yield a partner if the periodic orbit 
a; 7 lies close to a self-retracing path during one of the loops. This type of 'crossing 
region' is described by a 8y sr such that the original path x 1 and the time-reversed x l 7 
stay close together with \y^ r s \ <C 1 holding during the entire loop 1Z. The motion of the 
time-reversed path x^ in 1Z can then be obtained by linearisation around the original 
x 1 using the stability matrix R. In this case one finds R8y sr = F 5y sr neglecting 
corrections of higher then first order in y StU . The solution © is then replaced by 
8%R,i = 8y sr and Sx^ i = 0. It therefore does not give a new partner orbit but just 
the time- reversed periodic orbit x? = xL. But contributions to K(t) of this type 
are already treated in the diagonal approximation and must not be included in J2J). 
However, because of the hyperbolic nature of the dynamics the condition |y^ r J <C 1 
holds true for the entire loop 7Z only if the loop time T^ r is smaller than a certain 
minimal time 2fi im i n (x 7 (t), Sy), i.e. Tf[ < Tfj iin ; n . Corresponding arguments involving 
the time TL tia - ln (x 7 (t) , Sy) hold for the other loop C. These minimal times T^ ^) min 
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are determined by the time scale on which the linearization breaks down and are 
implicitly given by the condition 

y u ,s( ±T (R,L),min, X, Sy) = C UyS (x) (7) 

where the '+' ('— ') sign corresponds to y u (y s ) and c u , s (x) are constants of order one. 

Since © uniquely defines the partner orbit as sketched in figure |l(a)| one can 
show that the Maslov index for the orbit x 1 equals that of x?. The Maslov index of a 
periodic orbit is given by the winding number of the stable or unstable manifold ^21 • 
Since the partner orbit is close to the original orbit in section 1Z the contribution to the 
winding number accumulated between the two 'crossing regions' is the same for both. 
The second contribution comes from loop C and is thus given by the geometry xl^. 
The total winding number for x v , the sum of these two contributions, can be related to 
the total winding number of x in the following way. Time reversal symmetry implies 
the relation w s , u (x 1 (t)) — —Fw UyS (xl ( (T 1 — t)) between the manifolds of x 7 and x*. 
With this relation one can show that the contributions to the winding number coming 
from xly and x 7 during C are equal. Then the equality of the Maslov indices of x 1 and 
x? becomes evident |14) if one uses similar arguments as in the proof of the equality 
of the Maslov indices of a periodic orbit and its time reversed counterpart 

Action difference. The geometry of the partner given by JBJ allows one to derive 
the action difference between the two orbits of the pair (a; 7 ,cc^) as a function of Sy. 
Since the distance Sy is assumed to be small it is sufficient to expand the action in 
Sy. However, the expression one obtains for S ^H] still contains all the elements of 
the stability matrices R and L because the geometry of the partner orbit as given by 
(0 depends on them. Since the existence of the partner implies loop lengths larger 
than TrL,R),mm the vectors Sx have to lie very close to the respective local stable or 
unstable manifolds, e.g. Sxrj ~ y s w s (x). This fact enables us to express the action 
difference S^^ between the original orbit x 7 and the partner orbit x? = x^ in terms 
of the local manifolds [2] and the expansion coefficients y UiS given by (@J. Under 
the assumption that the directions of the manifolds are continuous functions of the 
position in phase space |§] the result then reads 

#7,7 = S(Sy,x) w (v^ Ziv s ) y s y u (8) 

which is correct up to second order in Sy. In the semiclassical limit it is sufficient to 
consider the regime given by |ytt,a| ~ y/h/i'Wu Zw s ) = y/h/S spo <C f which justifies 
the above restriction to small values of \y u ,s\- 

The equations JJjJ representing the geometry of x p and (JHJ are invariant under 
a shift of the PSS along the orbit within a 'crossing region'. This also implies that 
a 'crossing region' may include several selfcrossings in configuration space and hence 
the number of partners is not necessarily given by the number of selfcrossings as it 
was the case for the uniformly hyperbolic systems 

Counting 'crossing regions '. Moving the PSS along x 7 each 'crossing region' 
is characterized by a Sy that starts close to w s and ends almost parallel to w u . The 
number of 'crossing regions' can thus be determined by counting how often the unstable 
components y u of the vectors Sy go through a certain fixed value j/£ as one moves along 
Xj. This parameter y£ fixes the position within the 'crossing regions' used to identify 
and count it. The total number of 'crossing regions' must not depend on y£ and it 
will be shown that this is indeed the case. 
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To evaluate the weighted average (0 of the number Ng of contributing 'crossing 
regions' we make use of the following sum rule jll) valid for ergodic systems 

T-, T 

~Y,\^MT-Tj[dtf(x y (t))n Jdtf(x(t)) for T^oo.(9) 

7 

It relates the weighted average of a function f(x) over all periodic orbits of length T 
to a time average over a generic ergodic trajectory x(t, Xq) starting at any Xq in phase 
space. To apply © in the calculation of 10 one writes the number of events where 
the time reversed path comes close to the original one as 

d 2 iVg, 7 (a; 7 ft)) _ c ^ y u {0; Xl (t),S, y c u ) 

~~ d5d* wM*). »«V(<Wt),s,i/s)' (10) 

where p(x 7 (t); S, y°) is the density of partners per action in the PSS located at x 7 (t). 
The ratio y u /Vu describes the velocity of the flow in the PSS so that I|1U|I indeed gives 
the number of partners per action and time. In expression (|10H the position in the 
PSS is specified in terms of S and y° using ©. 

The density p(x;S,y^) is determined by the probability that the time reversed 
path goes through the point (S, in the PSS defined at x. The long time limit is 
thus given by the ergodic density p = T /12(E) = t/(2-kK) with 12(E) = (2irh) 2 d(E) 
being the volume of the energy surface in phase space. However, since certain 'crossing 
regions' with short loop lengths characterized by T(r^ l ^ < T( fl L ) !inin , see J7J), do not 
yield a partner one has to exclude parts of the time reversed path of length 

T min (a;, S) = 2[T R , min (x, S, y c u ) + T L , min (x, S, i/J] (11) 

with the factor 2 coming from time reversal symmetry. According to definition (JTJl 
this time T m - ln (x, S) is determined by the time it takes for the unstable component 
y u to grow from the small value S/(S spo c s ) <C 1 to the value c u ~ 1. It therefore 
does not depend on y£. To compute © we first apply © to i[TU)l. Then the minimal 
time l|llf) is of the same order as the Ehrenfest time and in the regime of small action 
differences S ~ h given by the asymptotic expression 

c u (x)c s (x)? s] 



,(x,S) w 2 In 



>1 (12) 



S 

where we used the standard definition of the Lyapunov exponent A in terms of any 
long ergodic path j^j . Using this relation we find that the main contribution to in 
the semiclassical limit is given by the averaged number of 'crossing regions' 



dS 



= pa 

7 



XT + 2 In const 



(13) 



where const depends on the structure of the phase space of the considered system but 
is independent of the action difference S and the time T The first of the three 
contributions on the r.h.s. of l|13|) is the largest one (~ h~ 2 ) and represents the ergodic 
properties of the system. The second term ~ h^ 1 In h reflects the underlying dynamics 
of the system [T^] and is much smaller than the first but still logarithmically larger 
than the third term ~ fir 1 . Inserting Ijl3(l into (J2J one finds that it is this logarithmic 
correction to the ergodic behaviour that gives the RMT result 

K^(t) w -2t 2 . (14) 

GOE-GUE transition. The crossover between the universality classes as time 
reversal symmetry is broken has been originally semiclassically obtained in Ref. |17| 



Leading off-diagonal contribution to the spectral form factor of chaotic quantum systems! 

within the diagonal approximation. Here we summarize a dynamical evaluation of this 
transition for the first off-diagonal correction. The appropriate transition parameter 
a interpolating between GOE and GUE is given by the ratio between the root mean 
square of a typical time-reversal symmetry breaking matrix element and the mean 
level spacing ^7]. We shall consider the case where the symmetry is broken by a 
uniform magnetic field B perpendicular to a uniformly hyperbolic two-dimensional 
system. Since for an orbit pair (x 7 ,x^) one of the common loops (£ or TV) between 
the 'crossing regions' are traversed in opposite direction in configuration space, the 
orbit pair acquires, owing to the overall magnetic flux enclosed, an additional action 
difference AnAB/fa. Here A is the enclosed (directed) area of the loop and 0o = h/2e 
the flux quantum. For hyperbolic systems the distribution of directed areas enclosed 
by trajectories of length t is, to good approximation 4 , Gaussian with variance t(3 
where /3 is a system specific parameter. To compute K^ 2 \t) at finite B we must 
additionally integrate, for given loop length t, the flux-induced action differences over 
the Gaussian area distribution. This results in a further damping exp[— t/tg] with 
l/t B = 2/3(2ttB/(/)o) 2 - Counting the ' crossing regions' with this additional weight 
eventually gives, together with the diagonal term, |14l I18j 

K(t; a) w t[1 + (1 - 2T) e - 8 ^ a2T ] for r -> (15) 

with a 2 r = (B I cf>o) 2 j3T . This precisely coincides with the form factor of parametric 
RMT [HI in the short time limit. 

To conclude, we have shown how correlations in the action of classical paths 
determine the spectral statistics of the quantum mechanical energy eigenvalues. We 
derived the next to leading order contribution for the spectral form factor K(t) 
and showed that it is identical to the corresponding RMT result. Our derivation 
is canonical invariant and based on phase space arguments only. This leads to the 
result that correlations in the classical action in hyperbolic chaotic systems are caused 
by 'crossing regions' in phase space where an orbit and its time-reversed version come 
close together. Since our method avoids the concept of crossings in configuration space 
it is suited to be extended to systems with more than two degrees of freedom. 
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